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Multi -Tethered  Space -Based  Interferometer:  Particle Jgyjatem  Model 


Introduction 

This  report  develops  a  dynamics  model  for  a  system  of  point  masses 
interconnected  by  idealized  tethers.  The  model  is  intended  for  the 
dynamic  analysis  of  certain  configurations  of  space  based 
interferometers.  The  interferometer  concepts  under  consideration  are 
comprised  of  multiple  collector  elements  and  auxiliary  bodies, 
interconnected  in  various  arrangements  by  low  mass  flexible  tethers.  The 
constellation  of  components  is  to  perform  coordinated  single  axis 
rotational  scanning  accompanied  by  baseline  dimensional  changes.  Spin 
axis  reorientation  and  spin-up/ spin- down  maneuvers  are  also  required. 
Figure-1  illustrates  several  candidate  interferometer  configurations 
which  fall  within  the  purview  of  the  model  developed  here. 

The  mechanical  system  considered  in  this  report  is  comprised  of  a 
system  of  point  masses  interconnected  in  an  arbitrary  fashion  by 
idealized  tethers.  The  tethers  are  treated  as  massless  tensile  members 
capable  of  exerting  force  only  along  the  straight-line  connecting  the 
respective  end  masses.  The  tethers  do  not  support  compression  or  any 
components  of  shear  force  or  bending  moment.  Both  extensible  and 
inextensible  tethers  are  considered.  The  constitutive  character  for  the 
extensible  tethers  is  taken  as  visco-elastic,  allowing  an  intrinsic 
energy  dissipation  mechanism.  The  amount  of  unstretched  tether  deployed 
between  end  masses  is  allowed  to  vary  in  a  specified  manner,  permitting 
study  of  deployment  and  retraction  as  well  as  fixed  length  operations. 
Each  particle  is  subject  to  an  arbitrary  force  arising  from  sources 
external  to  the  system.  Since  the  system  is  intended  to  operate  at  the 
L2  Lagrange  point,  gravitational  and  other  environmental  forces  are 
assumed  to  be  of  secondary  importance  relative  to  inertial  effects,  and 
hence  are  not  specifically  considered  here. 

The  dynamical  behavior  of  the  system  is  described  in  a  sequence  of 
analyses.  First,  the  general  nonlinear  unconstrained  motion  of  the 
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system  is  considered.  Important  system  linear  and  angular  momentum 
quantities  are  defined  and  Lagrange's  equations  are  used  to  produce  the 
motion  equations.  The  formulation  is  then  generalized  by  the  addition  of 
rheonomic  constraints  to  allow  for  specification  of  any  or  all  the 
degrees  of  freedom.  Lagrange's  equations  are  augmented  by  Lagrange 
multipliers  to  incorporate  the  constraints.  An  algorithm  is  provided  for 
the  solution  of  the  constrained  system  that  guarantees  satisfaction  of 
the  constraints  and  determination  of  the  associated  constraint  forces. 
The  unconstrained  nonlinear  motion  equations  are  then  linearized  by 
considering  small  amplitude  motions  about  a  nominal  steady  rotation.  The 
explicit  linearized  equations  are  obtained  and  properties  of  the 
coefficient  matrices  relevant  to  stability  analyses  are  established.  The 
nonlinear  dynamics  of  the  system  with  inextensible  tethers  and 
prescribed  motion  is  treated  in  the  Appendix.  Tether  inextensibilty  is 
introduced  through  the  addition  of  nonlinear  kinematic  constraint 
equations.  Equations  are  established  to  determine  the  tether  tension 
forces  and  any  additional  external  forces  needed  to  enforce  prescribed 
motion. 

Important  applications  of  the  models  developed  here  include: 
solution  of  particle  position  time  histories  relative  to  intended 
motions,  determination  of  control  forces  necessary  to  produce  intended 
motions,  and  evaluation  of  intrinsic  stability  characteristics  of  the 
system  under  conditions  of  steady  motion.  For  a  particle  system  free  in 
space,  which  is  intended  to  rotate  coherently,  and  for  which  analyses  of 
the  type  just  mentioned  are  of  interest,  the  choice  of  reference  frame 
in  which  to  describe  the  motion  requires  some  consideration.  Since  there 
is  no  rigid  body  in  which  to  embed  a  reference  frame  the  usual  "body 
fixed"  frame  of  attitude  dynamics  is  not  available.  Of  course  an 
inertial  frame  is  available,  and  is  indeed  necessary,  but  coordinates 
referred  to  it  are,  perhaps,  not  the  most  convenient  in  which  to 
linearize.  A  "mean  axes"  frame111*  could  be  defined  (provided  the  system 
possesses  a  sufficient  number  of  particles) ,  however,  its  utility  is 


*  Superscript  numbers  enclosed  in  brackets  indicate  reference  numbers. 
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diminished  if  small  vibration  modal  coordinates  are  not  introduced121 . 
Given  our  dynamical  interests,  the  application  of  a  mean  axes  frame 
appears  problematical.  These  considerations  have  lead  to  the  definition 
of  a  "prescribed  motion"  reference  frame  in  which  to  measure  the 
particle  positions.  The  rotation  of  this  frame  relative  to  inertial 
space  can  be  arbitrary,  but  is  assumed  to  be  fully  specified  as  a 
function  of  time.  It  will  be  found  that  the  use  of  this  prescribed 
motion  frame  facilitates  specification  of  prescribed  motion  studies  as 
well  as  linearization. 


Multi-Plane  Ring 


Figure-1.  Tethered  Interferometer  Concepts  within  Purview  of  Model. 
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System  Definition 


The  system  is  comprised  of  N  point  masses  PL  (i=l , 2 , . . . , N) .  In 
general,  Pi  can  be  connected  to  Pj  ,  (i*j),  allowing  the  possibility  of 
up  to  -^N(N-l)  connecting  tethers.  Let  be  an  inertial  reference 
frame,  with  origin  at  point-I,  and  unit  basis  vectors  I,  J,  K  .  5?p 

designates  a  "prescribed  motion"  reference  frame.  The  origin  of  ?p  is 
point-I,  and  its  unit  basis  vectors  are  denoted  ip,  jp,  kp .  The  rotational 
motion  of  ?p  with  respect  to  is  considered  to  be  a  known  function  of 
time.  Pi  has  mass  ir^  ,  and  position  vector  with  respect  to  point-I  given 
by  Pi  .  The  total  system  mass  is 

N 

=  X  mi 

i  =  l 

The  position  of  the  system  mass  center  with  respect  to  point-I  is 
denoted  p0  ,  and  is  defined  by 

1  N 

P®  =  -  Y,  miPi  (!) 

T~i 


Figure- 2  illustrates  the  system  geometry. 


.  4 


Figure-2.  System  Geometry. 


Let  v  denote  the  time  derivative  of  a  vector  v 
The  absolute  velocity  of  Pi  can  then  be  written  as 


observed  in  *PT  . 


Pi 


Pi  +  wxpi 


(2) 
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where  fi)  is  the  angular  velocity  of  5^  with  respect  to  5^  ,  and  pi 
designates  the  time  derivative  of  pA  observed  in  ?p  . 

Of  general  dynamical  interest  are  the  system  linear  momentum  and 
its  angular  momentum  about  point-I.  The  linear  momentum  of  the  system  is 
defined  as 

N  N  .  N  0 

L  =  ^  =  £  nUPi  =  ^  m±(  pA+  (0  X  p^  (3) 

i  =  l  i=l  i  — 1 

The  angular  momentum  of  the  system  about  point-I  is  defined  as 

N 

H=XPiXLi  <4> 

i=l 

If  the  vector  resultant  of  all  the  external  forces  acting  on  the 
particles  of  the  system  is  zero,  then  it  follows  from  Newton's  laws  that 

L  =  0  (5) 


Similarly,  if  the  vector  sum  of  the  moments  about  point-I  of  all  the 
external  forces  acting  on  the  system  is  zero,  then 


H  =  0 

The  kinetic  energy  for  the  system  is  defined  as 

N 

T  =  Pi 

i  =1 


(6) 


(7) 


Generalized  Coordinates 

The  system  configurations  of  interest  for  this  report  all  involve 
a  planar  arrangement  of  collector  elements  that  are  to  rotate  about  an 
axis  perpendicular  to  that  plane.  Radial  motions  of  the  elements  are  to 
accompany  the  rotations.  To  ease  the  prescription  of  such  motions  and  to 
facilitate  later  linearization,  cylindrical  coordinates  measured 
relative  to  ?p  are  selected  to  describe  the  positions  of  the  point 
masses.  The  cylindrical  coordinates  are  defined  such  that 
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Pi  =  r4  cos  0i  ip  +  rA  sin  0i  jp  +  zA  kp 


For  this  choice  of  coordinates,  it  should  be  evident  that  for  the  study 
of  the  nominal  rotational  motion  described  above,  kp  is  the  preferred 
axis  of  rotation.  Let  p.  denote  the  3x1  column  matrix  of  the  scalar 
components  of  pA  resolved  in  5?p  .  In  terms  of  the  vectrix  notation  of 


Hughes U1  we  have 


Pi  =  ?PT  =  ?PT 


f  rA  cos  0^ 
rA  sin  0A 

Zi  y 


(8) 


The  angular  velocity  (0  resolved  into  components  referred  to  ?p  ,  is 
expressed  as 


0) 


(0  = 


COy 


The  absolute  velocity  of  Px  given  by  Eg.  (2)  can  now  be  written  as 


Pi  =  <ii +  ®*Ei>  =  Z 


iL  cos  -  r^  sin  +  oo^  -  (0zr±  sin  0A 

ri  sin  0i  +  riQi  cos  0A  +  co^  cos  0A  -  (oxzx 
z±  +  co^  sin  di  -  (Oyri  cos 


(9) 


where  we  have  introduced  the  cross  product  operator 


b  = 

A  > 
bl 

=>  bx  = 

0  — b3  b2" 
b3  0  -bx 

b3 

k  3) 

-b2  bx  0  _ 

It  is  desirable  to  refer  certain  vector  quantities  to  their  scalar 
components  in  ^  .  Define  [CIp]  be  the  direction  cosine  matrix  that 
transforms  vector  components  referred  to  ?p  ,  to  those  referred  to  5^ . 

In  terms  of  the  notation  of  Hughes,  we  write 


=  rcIp]  ?, 


(10) 
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Let  the  vector  quantities:  p®  ,  ,  L  ,  and  H  have  scalar  components 

referred  to  5p ,  denoted  as  p0  ,  Li  ,  Lp  /  and  Hp  ,  then  we  can  write 


P®  =  ;  -i  =  mi<Ei  +  ^XEi) 

i  =  l 

N  N 

Lp  =  X  Hi  ;  Hp  =  X  Ei  Ei 


i  =  l 


i  =  l 


Resolving  vectors  p®  ,  L ,  and  H  into  components  in  7i .  and  denoting  the 
associated  3x1  column  matrices  by  R© »  Lx  /  and  Hj  respectively ,  we  have 

R®  =  t  CIp]  pffl 

hr  =  hP 
=  ICIp]  Hp 


The  quantities  Lj  and  Hx  are  conserved  under  the  conditions  described 
preceding  Eqs.(5)  and  (6). 

Substituting  Eq. (9)  into  (7)  yields  the  scalar  form  of  the  system 
kinetic  energy 

T  =  2  X  "it 

i  =  l 

+  Zx  [  zx  +  2 rA(  <ox  sin  0±  -  (0y  cos  0j  )  ] 

+  rf  [  <n?  +  ( odx  sin  0A  -  ffiy  cos  0x  f  ]  +  z?  (  ©x  +  <°y  ) 

-  2riz1toz(  <ax  cos  0J  +  (fly  sin  )  } 


r*  [  rt  -  2za(  cdx  sin  0x  -  CDy  cos  0x  )  ] 

+  r^Oi  [  r^  0A  +  2  <o2)  -  2  z^  (0X  cos  0A  +  (Oy  sin  0i  )  ] 


Lagrange  *  s  Equations 

The  equations  of  motion  are  first  developed  assuming  the  3N 
generalized  coordinates  are  independent.  Later,  prescribed  motion 
constraints  are  introduced. 

Assuming  the  generalized  coordinates  ri#  0i#  zi  ,  (i=l ,  2  ,  .  .  . ,  N)  are 
independent,  Lagrange's  equations  for  the  system  can  be  written  as 
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(12a) 


d 

dt 


3t 

3rA 


Qri 


d  (jto  "I  dr_  = 

dt  aei  90i  Qei 

V  1  / 


d 

dt 


3t 

aZi 


Qzi 


(12b) 


(12c) 


for  i=l, 2, .  . . ,N.  Qri  ,  Qei  «  and  Qxi  are  the  generalized  forces  associated 
with  the  indicated  degrees  of  freedom,  arising  from  both  internal  and 
external  forces  acting  on  the  particles.  Performing  the  indicated 
derivatives  on  the  left  hand  side  of  Egs. (12)  we  record  the  following 


d  (  3t  > 

dt  dr i 
V  1  / 


Hi*!  -  mi(£J2iZi+  ^2iZi> 


(13a) 


d  fdr  ' 

dt^e^ 


*irf 0A  +  mi[ri©z+  2riri(6i  +  coz)  -  £2li(rizi+  riZi)  -  ^£2^] 


(13b) 


d  fdT  ' 

dt  dzt 
V  1  / 


triiZi  +  mi(£22ir1 +£2^^) 


(13c) 


—  =  tni[(ri0i  +  2ri(0z-  £2^)  ^  +  £2*^  +  (  toz  +  £2^  )  rA  -  coz£2lizi ]  (14a) 

— -  =  rn. A [  £2li(r^zi  —  ri  ^2i )  ©z^  ^  (14b) 

O 

=  mj  (©*  +  C0y)  Zi  -  0i  +  <0Z)  “  ^2i^i  ]  (14c) 

3z± 

In  Eqs. (13)  and  (14)  we  have  introduced  the  following  quantities: 


=  0)x  cos  0i  +  (Oy  sin  0i  (15a) 

Cl2i  =  G)x  sin  0A  -  0)y  cos  0i  (15b) 

an  =  6)x  cos  0A  +  d)y  sin  0A  (15c) 

«2i  =  ©x  sin  “  ©y  COS  0i  (15d) 

=  axi  -Q2i0i  (15e) 
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<15f) 


^2i  —  a2i  +  ^li®i 

Generalized  Forces 

The  forces  acting  on  the  particles  may  be  divided  into  those 
arising  from  sources  external  to  the  system  and  those  arising  from  the 
interconnections  between  particles.  Consider  first  the  external  forces 
acting  on  P±  .  Let  be  the  resultant  of  the  external  forces  acting  on 
p  .  it  proves  convenient  to  resolve  FA  into  three  orthogonal  components 
along  directions  that  conform  naturally  to  the  cylindrical  coordinates. 
For  Pi  we  define  a  Cartesian  reference  frame  ,  having  unit  basis 
vectors  e^,  e(ti},  where;  is  directed  along  (Pi“Pi*kp),  e^  is 

parallel  to  kp  ,  and  e^  =  e™  x  eJH  .  Then  we  can  write 

Fri 

Fti 
FZi 

Figure- 3  illustrates  the  situation  in  the  ipjp  plane. 

7>i 

y\ 

i 

Figure- 3.  Resolution  of  External  Force  on  Pi  . 

From  Eqs . (8)  and  (10)  we  have 

Pi  =  5?  lClp(t)  ]  p.(ri(  0i(  zA ) 

The  virtual  displacement  of  is  given  by 
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where  8^ ,  80A  ,  and  bzi  are  the  virtual  displacements  of  the  generalized 

coordinates  of  Pi  .  The  direction  cosine  matrix  transforming  vector 
components  from  5?  to  is  given  by- 

cos  sin  0A  0 

[  cipl  =  -sin  0A  cos  Qi  0 

0  0  1 

Resolved  into  scalar  components  referred  to  5^  ,  8pi  can  be  expressed  as 

r  8*i  ' 

8Pi  =  9i 

<  8zi  > 

The  virtual  work  performed  by  the  external  forces  acting  on  the  system 
is  given  by 

8w(e>  =  £  Fi  •  8p± 

i=l 

N 

=  X  (  Fri8ri  +  Ftiri60i  +  Fzi8zi) 
i=l 

=  X  (  Q S  fci  +  Q'ei’  89 1  +  qS  8Zi )  (18 ) 

i=l 

In  Eq.  (18)  we  have  introduced  the  notation,  ,  to  designate  the 
generalized  force  associated  with  generalized  coordinate  , 
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(Pi  =  r-i.e^Zj)  arising  from  the  action  of  the  external  forces  acting  on 

Let  f1;i  be  the  force  acting  on  Pi  arising  from  its  connection  to 
p  To  allow  arbitrary  connectivity  we  assume  that  each  particle 

could  be  connected  to  every  other  particle  of  the  system.  The  vector 
resultant  of  all  internal  forces  acting  on  Pi  is 

fi  =  Xfid 

j=l 


The  virtual  work  performed  by  f*  on  Pi  is  given  by 

8wf ’  =  fi  •  Spi  =  2  fid  • 
d=i 


The  virtual  work  performed  by  all  the  internal  forces  acting  on  all  the 
particles  of  the  system  is 


8wa>  =  X  X  ?id  '  ^i 


i— 1 j=l 


(19) 


By  assumption  the  tethers  are  only  capable  of  exerting  tensile 
force  along  the  straight  line  connecting  their  respective  end  masses, 
thus 


fid  =  fid  £  0 


Let  li:j  be  the  position  vector  from  Pt  to  Pj ,  i.e. 


1ij  “  Pd  “Pi  “ 


Xj  COS  0J  -  rA  COS 

rj  sin  0^  -  r*  sin  Bj 


Zj  -  zA 


The  distance  between  Pi  and  Pi  is 


^id  “  Iji  ”  lij 


(20) 


(21) 


+  rj1  -  2riacj  cos  (Bj  -  Bi)  +  (  Zj  —  z*)2 


(22) 


11 


and  its  time  derivative  is 


iij  =  I  ri  ii  +  K)ij  -  (*i*j  +  riEj)cos(0j-ei) 

+  riEj(  Qj  -6i)  sin  (  Bj-Oi)  +(  Zj  -z*)  (  Zj  -z^  ]  /  ly 

To  allow  the  possibility  of  energy  dissipation  in  the  tethers, 
their  constitutive  character  is  assumed  to  be  visco-elastic.  Further,  we 
assume  that  the  length  of  unstrained  tether  deployed  between  any  two 
particles  of  the  system  is  a  known  function  of  time.  Let  the  scalar 
dij(t)  be  the  length  of  unstrained  tether  deployed  between  Pi  and  P^  at 
time  t.  Temporarily  dropping  the  particle  subscript  notation,  the 
longitudinal  strain  of  a  tether  is  defined  as 

l(t)  -  d(t) 

e  =  - 

d(t) 

Assuming  a  Kelvin- Voigt 143  constitutive  law,  the  tensile  stress,  T  ,  and 
the  longitudinal  strain  are  related  by 

X  =  e(e  +  ae) 

where  E  is  the  effective  Young's  modulus  for  the  tether  and  a  is  a 
constant  dissipation  parameter.  If  A  is  the  effective  cross  sectional 
area  of  the  tether,  then  the  tension  is  given  by 

\  AX  £  >  0,  T  >  0 

H  o 


Assuming  that  e<<l,  and  re-establishing  the  particle  subscript 
notation,  the  tension  in  the  tether  connecting  Pi  and  P ’■  can  be  written 
as 


=  < 


kijdij-dij)  +  Cydij-dy)  if 


otherwise 


li^dy 

cij(dij  —  di;j)  ci  — ki;j(li:j—  di;j) 


(23) 


where  the  effective  stiffness  ky  ,  and  damping  coefficient  ci;j ,  are 
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kij  -  kji 


cij  -  cji  - 


(  EA  ) 


ij 


*ij 


The  virtual  work  done  by  the  internal  forces  of  the  system,  given 
by  Eq.  (19),  can  now  be  expressed  as 


i=l j=i 


Using  Eqs.(17)  and  (21)  we  get 

8w(i)  =  £  [  Qg  Sr,  +  80,  +  Q(Ji  8z,  ] 

i=l 

where  the  generalized  forces  associated  with  the  internal  forces  are 
given  by 


II 

§i  = 

Qei 

Qzl 

\  ) 

|  5-i  1 

,  j 

In  Eq. (24)  we  have  introduced  the  quantities 

C- j  =  3Cj  cos  ( 6j  -  0A)  -  r± 

'Cij  =  -"Cji  =  ZiZj  sin  ( 6-j  6, ) 

Cij  =  "Cji  =  <=S-Zi) 


(24) 


(25a) 

(25b) 

(25c) 


The  generalized  forces  appearing  on  the  right  hand  side  of  Eg. (12)  can 
now  be  written  as 

Qpi  =  Qpi  +  Qpi  (26) 

(pi  =  ri#  0i,  zA)  . 

Substituting  Eqs.  (13),  (14)  and  (26)  into  Eqs.  (12)  yields  the 
nonlinear  motion  equations  for  the  system 

[M],u,  =  Q(,e)  +  N,  +  S,  (27) 


13 


(  i  =  1,  2,  .  .  . ,  N  ) ,  where 


’  1 

0 

0  ’ 

(  V  \ 

*ri 

[Mli  =  mi 

0 

ri 

0 

'•  Hi  = 

01 

;  Q(.e)  = 

riFti 

0 

0 

1 

F  • 

1  21  ; 

— i 


)  +  z^aji-to^a)  +  ^(^  +  20^)  ^  +  2Q2iZi 

rf(  ila^zi  -®z  )  +  r^J  Qu  +  ^2i<  ®i  +  )]  +  2x^ -  ( ^  +  (0, )  !•  ] 

z^  a£  +  oo£ )  -  2£22i  r-  -  rA  t  ctji  +  Qu(  2  Qx  +  coz )  ] 


(28) 


Prescribed  Motion 


Define  the  vectors  of  generalized  coordinates  and  forces  for  the 
system  as 

{  <1  }T  =  {  ri  ®l  Z1  r2  Z2  *  •  *  rN  ZN  }  =  {  ^2  *  *  *  ^3N  } 

{  Q  }T  =  {  Qrl  Qei  Qzl  Qr2  @02  @z2  *  *  *  @rN  @GN  @zN  }  =  {  @1  @2  •  *  *  @3N  } 


Using  the  above  notation  and  that  of  Eq. (26) ,  Lagrange's  equations  for 
the  unconstrained  system  can  be  written  as 


+  Q 


(i) 

k 


for  k  =  1,  2,  .  .  . ,  3N  ,  or,  collectively,  as 


d 

dt 


(29) 


For  the  study  of  prescribed  motion  we  assume  that  np  of  the 
generalized  coordinates  {q},  are  known  functions  of  time.  We  assume  the 
presence  of  rheonomic  constraints  of  the  form 

{qp(t)}  =  [Sp]T{q}  (30) 

where  {qp(t)}  is  an  npXl  column  matrix  of  known  time  functions  (  np  <  3N  )  , 
and  [ Sp ]  is  a  3NXnp  Boolean  selection  matrix.  The  generalized 
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coordinates  have  been  chosen  so  that  the  system  motions  of  interest 
could  be  described  in  this  simple  manner.  It  is  known  from  analytical 
mechanics  that  Lagrange's  equations  for  a  system  subject  to  constraints 
of  the  form  given  by  Eq.  (30),  can  be  written  as 

where  {X}  is  an  npxl  vector  of  Lagrange  multipliers. 

Using  the  nomenclature  established  for  Eq. (27) ,  the  nonlinear 
motion  equations  for  the  system  subject  to  the  constraints  of  Eq. (30) 
are  now  written  as 

lW]{q)  =  {Q(e)}+{«}  +  {S}-[Sp]{M  (32) 

where 


'  [Mh 

[0] 

...  [0] 

Q(e)  "l 
-1 

t  an  = 

[0] 

[M]2 

...  [0] 

{ Q<el  }  = 

Q(e) 

—2 

_  [0] 

[0] 

-  [M]n 

Q(e» 

l  -N  J 

'  §i  ' 

m  = 

n2 

t 

{5}  = 

§2 

Nk 

K  -*  ) 

The  last  term  on  the  right  hand  side  of  Eq.  (32)  can  be  interpreted 
as  the  generalized  force  required  to  enforce  the  constraints,  .i.e. 

{ Qfc>  }  =  ~[Sp]{X}  (33) 


To  establish  the  physical  forces  of  constraint  some  specific 
deterministic  arrangement  of  scalar  force  components  must  be  considered. 
Here  we  shall  assume  a  system  of  constraint  forces 

Ff}  ,  i  =  1,  2,  .  .  . ,  N  ,  (not  necessarily  all  nonzero)  which  act  in  parallel 
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with  the  applied  external  forces  defined  by  Eq. (16) .  Resolved  in  frame 


5^  the  components  of  Ff0  are  designated  as 


5(0  _ 


= 


/  T7(c)  N 
*ri 

T?(c) 

p(c) 

v  2i  J 


Equation  (18)  specifies  the  relations  between  the  components  of  the 
external  forces  and  the  corresponding  generalized  forces.  Relations 
identical  to  those  of  Eq.  (18)  exist  between  the  components  of  FiC)  and 
the  generalized  constraint  forces  embodied  by  Eq.  (33)  .  Thus,  the 
physical  constraint  forces  are  related  to  the  Lagrange  multipliers  by 


=  -rsPHM 


(34) 


Equations  (3  0)  and  (32)  constitute  a  system  of  3N  +  np 
differential-algebraic  equations  in  3N  +  np  unknowns,  { q }  and  {X-}.  Given 
the  simple  form  of  the  constraint  equations,  Eq. (30) ,  a  coordinate 
transformation  can  be  made  which  yields  a  set  of  differential  equations 
in  terms  of  an  independent  set  of  generalized  coordinates,  which  are 
uncoupled  from  the  Lagrange  multipliers.  Define  the  set  of  generalized 
coordinates  which  are  not  prescribed  (i.e.  which  are  "free”)  as 

{q£}  =  tSf]T{q}  (35) 


where  [  Sf  ]  is  a  3Nx(3N-np  )  Boolean  selection  matrix.  Equations  (30)  and 
(35)  can  be  written  together  as 


r  {qf}  1 
{qP} 


t  S]T{q  } 
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It  follows  from  the  simple  Boolean  structure  of  the  matrices  [Sf]  and 
[Sp],  and  the  compl ementary  nature  of  the  coordinates  {qf}  and  {qp}, 
that  [S]  is  an  orthogonal  matrix.  Consequently,  we  have 


{q}  =  [rs£]  [Sp]] 


{qP} 


(36) 


[SD]T  [SD]  =  [1]  npXnp 


[Sf]T  [SD]  =  [0]  (3N-np  )xn„ 


Substituting  Eq. (36)  into  Eg. (32),  and  premultiplying  the  result  by 
[S]T,  yields  the  partitioned  system 

[Sf]T[WHSf]{qf}  =  [Sf]T  ({Q(e)} +{*}  +  {*} -[*]  [Sp]{qp})  (37) 


{Jl}  =  [Sp]T({Q(e,}  +  {«>  +  {^}-t^](tSf]{q£}  +  [Sp]{qP}))  (38) 


Given  {qp} ,  {qp},  and  {qp},  Eq.  (37)  can  be  integrated  to  determine  {qf}, 
and  Eq.  (38)  can  be  evaluated  to  determine  {A^t)}. 


Linearization  of  the  Unconstrained  Motion Equations 

In  this  section  we  linearize  the  unconstrained  motion  equations 
about  a  state  of  steady  spin  with  the  tethers  of  fixed  unstrained 


length.  Throughout  this  section 

0)  =  0  (39a) 

djj  =  0  i,  j  =  1,  2,  .  .  N  (39b) 

The  coordinates  of  JDj  relative  to  ?p  are  taken  to  be 

rA  =  rf  +  ArA(t)  (40a) 

0A  =  0j  +  AG^t)  (40b) 

zi  =  z\  +  AZ^t)  (40c) 
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where  r°,  0°,  and  zj  are  constants  defining  the  steady  state  equilibrium 
configuration.  The  variables  Ar^  A0if  Aza  and  their  time  derivatives  are 
considered  to  be  quantities  of  small  magnitude,  and  we  will  ignore  terms 
in  which  they  appear  to  second  or  higher  order.  Following  the  notation 
established  in  Eqs.(40),  we  introduce  also 

ui  =  u°  +  Au,  (t) 

{q}  =  {q0}  +  (Aq(t)} 


(41a) 

(41b) 
(41c) 

(41d) 

(41e) 

Above  we  have  introduced  the  constants 

=  ©x  cos0°  +  (0y  sin0° 

‘  ‘  =  cox  sin  0°  -  coy  cos0° 

Substituting  Eqs.(40)  and  (41)  into  Eq. (28)  and  ignoring  second  and 
higher  order  terms,  we  get 


Recognizing  Eq. (39a)  and  substituting  Eq. (40b)  into  Eqs. (15) 

yields 

=  Q°u  -  £2*1  A0i 
Q2i  =  +  £&  A0, 

“li  =  «2i  =  0 

£2U  =  -£2°x  A0j 

£^2i  — 


Nj  =  -[G]i  Au.  -  [C]i  Au.  +  Nj  (42) 

where 


'  Ar^ 

0 

cT 

A0i 

i 

[Glx  =  2mj 

<*>z*i 

0 

AZi 

V  1  ) 

«Si 

O0  T*0 
1Zli  ri 

0 
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t  C3± 


-[a^+(£4f] 


itii 


-Q^Q^rf+m.z?) 


-^[(c^+ffl^rf+ax^iZ?]  -c^Q^r? 

-to^r®  -«o£  +  <h£) 


m 


'  [^  +  (£2L)2]r®  -  to^uz?  ' 
£2°i(Q^rf+©zz0i)r10 
-®z  ^li  ri  +  (®i  +®y  )  z? 


The  linearized  counterpart  to  the  mass  matrix  appearing  in  Eq. (27)  is 


[Molt 


10  0 
0  (r?f  0 
0  0  1 


(43) 


Linearizing  Eqs. (25)  about  the  steady  state  solution  we  have 

®ij  =  <*?j  +  <*ij  <44a> 

Xij  =  t?j  +  tij  <44b> 

;«  =  ti»  +  tii  <44o 

where  the  constant  terms  are  defined  as 

<r?j  =  r?  cos  ( 0?  —  0?  )  —  rf  (45a) 

xlj  =  -x°n  =  sin (05-0?)  (45b) 

tij  =~tii  =  ~  *°i  <450 

and  the  terms  dependent  to  first  order  on  the  deflections  from  the 
steady  state  are  recorded  as 

O^j  =  cos  (0°  -0°)  Arj  -  Ar.  -  r?  sin  ( 0°  -  0?)  ( A0..  -  A0.)  (46a) 

tlj  =  rfsintOj-OjiArj  +  if  sin  ( 0?  -  0? )  Ar*  +  rf  Ej° cos  (0®  -  0?  )  (A0j  -A0J  (46b) 

=  AZj  -  Aza  (46c) 
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Expanding  Eq. (22)  to  first  order  we  have 


1  —  1° 
±ij  “ 


+  1id 


where 


and 


lij  =  4-  l-CTdi  Ard  -  Ari  +  T°i  (A0j-A0i)  +  Cij  (Az^-AZi)  ]  (47) 

1ij 


We  also  record 


(l..)"1  «  J: - 

1]  lij  dij  )2 


(48) 


and 


Lij 


=  iij 


=  4“  [-cji  A^j  -  aij  ArA  +  T°j  (AOj-AOi)  +  £y  (A^-Az^  ]  (49) 


For  our  stipulated  conditions,  Eq. (23)  simplifies  to 


fij  -  ^ij  +  ^ij  +  Cij 


where 


f°i  =  ky  (  1^  -  dy  ) 


(50) 


(51) 


is  the  steady  state  tension  in  the  tether  connecting  and  Pj  . 

Substituting  Eqs. (44) , (48)  and  (50)  into  Eq. (24) ,  and  dropping 
terms  above  first  order,  yields 


(52) 
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+ 


N 


I 


\ 


ii: 


/ 


Introducing  Eqs . (46) , (47)  and  (49)  into  Eq. (52),  we  can  write 

S±  =  S °  -[«r]i{Ag}-[r]i{Aq  } 


where 


'a° 


H  xrO 

s-x£ 

j=l  Xij 


ij 


[*rii  =  -  [  [k]a1  [k]j 


K]«  ] 


Introducing 


Krr 

Kre 

Krz 

[Wij  = 

Kflr 

Kee 

Kez 

Kzr 

KZ0 

Kzz_ 

ij 

/ 

fo  A 

0  _  0  _ 
eij  eji 

<*?j 

ra 

kij  " 

£jj 

3° 

Xij  J 

we  have  for  i  *  j  : 


(Krr  >ij  =  <Krr  Jji  =  ^  COS  (8!?  -  6?) 


Lij 


(  Krf  )y  = 


-0  «0  T0 

eij  °ij  X 


ij 


--f  tfBin(e;-e;) 

1ij 


(Krz  >„  = 

(  Kar  )ij  =  (  Krt  >ji  =  -e?j  Oji  x°ij  +  751  ri  8in  ( ej  -  0i) 

Xij 

(  Kee  Uj  =  (  Kee  )ji  =  e?j  ^Sj  ?  +  To^  ri  rj°  cos  (®j -  ®i) 

iij 


(53) 


(54) 
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(  K0z  )ij  ~  (  K0z  )ji  ~  eij  Tij  ^ij 
(  ^zr  ^ij  =  (  ^rz  ^ j  i  “  —  eij  ^ji  ^ij 

(  Kz0  )ij  =  (  Kz0  )ji  =  (  K0Z  )ij 


(Kz2  hi  =  (Kzz  )j4  =  ejj  (Cij  )2  +  -y- 


For  j  =  i  we  have: 


f° 

e“j(o?j)2+41  (Rre  >ij  <Krz)ij 
1ij  , 


(Kee)ij  <Kez)ij 


(Kzz  )y 


For  the  coefficient  matrix  of  the  last  term  on  the  right  hand  side  of 
Eq. (53)  we  have 

[Z>]i  =  -[[D]u  [D]^  -  [D]iN  ]  (! 


V>hi  = 

Uij  ) 


_0  _0 
-CTijCTji 

CTij  Tij 

_o  ro 

_0 

_Tij  °ji 

(xj,  f 

t°  r° 
Tij  *ij 

r  o  „.o 

(C^)2 

i  *  j 


Evidently  [D]^  is  not  symmetric.  We  have  also 


[Dili  = 


_  Cij 

A  (10  v2 
j=l  Uij  ' 


7^0  \2  _0  _0  _0 
<aij  )  °ijTij  aij  ^>ij 

<^>2 

sym  (CSj  )2 


Substituting  Eqs. (42) , (43)  and  (53)  into  Eq.(27),  and  rearranging, 


yields 


[M0  AUj  +  [G]i  Au.  +  [Cli  Au.  +  t  Z?  ]±  {  Aq  }  +  [^{Aq} 
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of  +  N°  +  S° 


for  i  =  1,  2,  .  .  . ,  N  .  Assembling  the  system  of  equations  given  by  Eq.  (56)  as 
a  single  matrix  equation,  yields 

[W0  ]  {Aq}  +  ([?]  +  IV)  ){Aq}  +  (  [  <?  ]  +  [  2T  ]  ){Aq} 

(57) 

=  {q<c1  }  +  {%}+  {^0} 

where 


'  [Moll 

[0] 

...  [0] 

'[G]x 

[0] 

...  [0] 

IW0]  = 

[0] 

[M0]2 

...  [0] 

;  [  ?  1  = 

[0] 

[G]2 

...  [0] 

[0] 

[0] 

-  [MoIn 

[0] 

[0] 

-  IG]n 

[Cli 

[0] 

...  [0] 

to] 

[C]2 

...  [0] 

[0] 

[0] 

-  [C], 

'  [t>\  ' 

'  [XTli  ' 

[Z>]  = 

[2>]2 

;  t  X  ]  = 

[*ri2 

.  . 

- 1 

S3 

»—i 

u 

_ 1 

ft  > 

r§n 

£ 

U }  - 

\  ”  ) 

S° 

^  “N  / 

The  vector  of  external  forces  appearing  in  Eq.  (57)  can  be 
separated  into  two  terms,  one  constant ,{  o[je)  }  ,  and  the  other  time 
dependent ,  {  Q^t) }  ,  i .  e . 

{  Q<e>  }  =  {  Qoe>  }  +  {  Qie,(t)  }  (58) 

If  we  consider  Eq. (57)  under  the  conditions 
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{  Aq  }  =  { Aq }  =  {  Aq  }  =  {  Q^’jt) }  =  {  0  } 


then  we  find  the  steady  state  equilibrium  equations  for  the  system  are 

{QS>e,}  +  {3U  +  Uo}  =  {0}  (59) 

Equations  (59)  constitute  a  system  of  3N  nonlinear  algebraic  equations 
in  the  quantities;  r ? ,  0°,  z°  ,  i  =  1,  2,  .  .  .  ,  N  . 

Recognizing  Eq. (59) ,  the  linearized  motion  equations  for  the 
unconstrained  system  follow  from  Eq. (57) ,  as 

[7X0]{Aq]  +  ([?]  +  IV]  ){Ag}  +  ([<*]  +  [5?]  ){Aq}  =  fae,(t)}  (60) 

The  following  properties  of  the  coefficient  matrices  appearing  in 
Eq. (60)  are  noted: 

m0]  =  m0f  =  [V]  =  [Vf 

i  e  i  =  t  e  ]T  [  x  ]  =  t  f 


Model  Sincmlaritv 

A  unique  solution  to  any  of  the  various  sets  of  motion  equations 
derived  in  this  report,  (e.g.  Eqs.(27),  (37),  or  (60)),  requires  that 

the  coefficient  mass  matrix  of  the  generalized  accelerations  be 
nonsingular.  For  the  unconstrained  systems  (linear  and  nonlinear),  the 
mass  matrices  [  ]  /  and  [  ]  ,  are  diagonal  and  possess  elements  rA  and 

xf  ,  (i=l ,  2 ,  .  . .  ,N)  ,  respectively.  Evidently,  if  ^=0  or  r°  =0  then  the 
associated  mass  matrix  would  be  singular.  This  is  a  consequence  of  our 
choice  of  coordinates  rather  than  an  indication  of  a  physical 
degeneracy.  This  deficiency  of  the  models  can  be  avoided  by  limiting 
attention  to  only  those  conditions  for  which  rA  *0  (or  r°  *0  for  linear 
analyses),  (i=l , 2 , . . . ,N) .  The  situation  with  the  constrained  nonlinear 
equations  is  somewhat  more  involved.  The  generalized  mass  matrix,  from 
Eq.  (37)  ,  is  [Sf]T[3£]  [Sf]  ,  and  it  can  be  shown  that  that  matrix  is 
nonsingular  if  [  W  ]  is  nonsingular.  ([Sf]T[flfc]  [Sf]  may  also  be 
nonsingular  when  [  Wt  ]  is  singular,  depending  on  which  degrees  of 
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freedom  are  specified.)  The  restriction  to  conditions  for  which  1^*0, 
however,  guarantees  a  unique  solution  for  Eq. (37)  also. 


Kinematics  of  7r 

The  motion  ?p  relative  to  has  been  assumed  to  be  a  known 
function  of  time.  To  complete  the  formulation  we  shall  define  the 
orientation  of  ?p  with  respect  to  %  in  terms  of  a  general  1-2-3  Euler 
angle  sequence.  The  orthogonal  unit  basis  vectors  i,  j,  k  for  each  frame 
referred  to  below  will  be  denoted  with  appropriate  subscript  and 
superscript  identifiers.  The  signs  of  all  single  axis  rotations  are 
positive  according  to  the  right  hand  rule.  Let  5^  be  an  intermediate 
frame  achieved  from  by  a  rotation  of  angle  about  the  I  axis.  Let 
frame  ?p  be  achieved  from  by  rotation  of  angle  i|/2  about  the  axis  jp  . 
?p  is  achieved  from  ?p  by  a  rotation  of  angle  ^3  about  the  kp  axis.  The 
direction  cosine  matrix  transforming  vector  components  from  ?p  to  ?x  is 
obtained  by  concatenating  the  sequence  of  single  axis  rotations  yielding 


ci|f2c\|r3  -C^2S^3  S^2 

cil/1BTjf3  +  s^si^cta  cn^ci^  -  si^s^sts  “S^c^ 
si|f1s\lf3  -  c^si^ci^  ST)f1cn|;3  +  c^s^si^  civets 


where  we  "have  used  the  shorthand  notation:  ci)^  —cos  ,  si)^  —  sin  i(ri  . 
The  angular  velocity  of  ?p  with  respect  to  5^ ,  expressed  in  terms  of 
the  Euler  angles  is 

5  =  I  +  iir2  jp  +  kp 
Resolved  into  conq?onents  referred  to  ?p  we  have 

*2 
4*3  _ 

=  [n]  \jr 


(0  = 


c4r2ci)/3  si)f3  0 
-ct2s4f3  ci|f3  0 
sijr2  0  1 
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Appendix  -  Inextensible  Tethers 


In  the  main  body  of  this  report  we  have  developed  the  equations 
governing  the  motion  of  a  system  of  particles  interconnected  by 
idealized  extensible  tethers  subject  to  applied  external  forces  as  well 
as  prescribed  motion.  The  formulation  presumes  that  the  length  of 
unstrained  tether  deployed  at  any  instant  is  known  a-priori.  Clearly, 
the  unstrained  length  of  tether  deployed  has  a  strong  influence  on  the 
tether  tension,  which  is  to  be  a  principal  means  of  stabilization  and 
control  of  the  system.  In  certain  circumstances  (e.g.  inverse  dynamic 
analysis)  the  need  to  specify  tether  unstrained  length  time  history 
proves  to  be  problematical  and  simplification  is  obtained  by  treating 
the  tethers  as  inextensible.  Inextensible  tethers  are  introduced  through 
the  definition  of  inter-particle  distance  constraints.  Unlike  extensible 
tethers,  for  which  tension  is  described  by  constitutive  equations  (as  in 
Eq.  (23) ) ,  with  inextensible  tethers  the  tensions  must  be  solved  for  as 
constraint  forces.  In  this  appendix  all  tethers  of  the  system  are  taken 
to  be  inextensible.  Supplemental  prescribed  motion  constraints  are  also 
considered  and  the  forces  necessary  to  produce  those  motions  are 
determined. 

Let  nt  be  the  total  number  of  tethers  in  the  system,  and  let 
tether  k  connect  particles  Pi  and  P^  .  The  distance  between  Pi  and  P^  is 
given  in  Eq.  (22)  by  li;j  .  Let  di;j(t)  be  the  length  of  inextensible  tether 
instantaneously  deployed  between  Pi  and  P ^  .  We  shall  assume  that  dAj(t) 
is  a  known  function  of  time.  If  Tk  is  the  tension  in  tether  k,  then  the 
constraint  governing  tether  inextensibility  is 


iij  =  dij(fc>  if  Tk  >  ° 
no  constraint  if  Tk  <  0 


(61) 


Whether  the  kinematic  constraint  of  Eq. (61)  is  in  effect  or  not  depends 
on  the  tension  force,  so  the  system  is  seen  to  possess  a  discontinuous 
number  of  degrees  of  freedom.  Rather  than  deal  with  the  complication  of 
the  discontinuity,  we  shall  simply  enforce  the  constraint  at  all  times 
and  recognize  that  the  occurrence  of  non-positive  tension  invalidates 
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our  results  from  that  point  onward.  We  note  that  taking  Eq.  (61)  to  hold 
for  all  Tk  corresponds  to  the  condition  of  a  rigid  link. 

It  proves  convenient  to  state  the  tether  constraints  in  the  form 

hj^r^  0i#  zif  rj,  0j,  Zj,  t)  —  { 1± j  )2  -  (d^  )2  =  0 
or 

hk(q,  t)  =  r42  +  ij2  -  ^Ej  costej-ei)  +  (z^-z  -  (di;j)2  =  0  (62) 


For  the  entire  system  of  nt  tethers  the  constraints  are  written  as 


{  h(q,  t)  } 


h-^q.t) 


h2(q,t) 


=  0 


hnt(q,t) 


(63) 


To  permit  the  study  of  prescribed  motion  and  to  allow  the 
determination  of  the  associated  control  forces  we  consider  in  addition 
to  Eq.  (63)  the  np  constraints  given  by  Eq.  (30)  .  The  complete  set  of 
nt  +np  constraints  is  expressed  as 


{<I>(q,t)  } 


{  h(q,t)  } 

[Sp]T{q}  -  {qp(t)} 


=  0 


(64) 


It  is  required  that  nt+np  £  3N ,  and  that  the  constraints  of  Eq.  (64)  be 
consistent  and  functionally  independent. 

Lagranges  equations  for  our  system  of  particles  subject  to  the 
constraints  of  Eq. (64)  can  be  written  as 


which  must  be  considered  together  with  the  constraint  equations  given  by 
Eq.  (64) .  Note  that  in  comparison  with  Eq.  (31) ,  Eq.  (65)  does  not  contain 
the  generalized  forces  associated  with  any  internal  forces,  since  here 
the  tether  tensions,  being  constraint  forces,  do  no  work  in  virtual 
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displacements  consistent  with  the  constraints.  In  Eq. (65)  we  have 
introduced  the  Jacobian  of  the  constraints 


where 


]  = 

r— i 
W 

[Spf 

- 

dX 

3hx 

__ 

3hj 

_ 

3hi 

ar2 

3hx 

3zn 

'  aiT 

3h2 

3h2 

30! 

3h2 

3z! 

3h2 

3^7 

3h2 

dzN 

,3q. 

3hnt 

3hnt 

^nt 

*lnt 

_ 

30! 

3z! 

3r2 

3zn 

Forming  the  partial  derivatives  indicated  above  from  Eq.  (62)  we  record 


for  l  =  8i:i  is  the  usual  Kronecker  delta  symbol.  In  terms  of 

the  notation  introduced  with  Eq. (32),  Eq. (65)  can  be  written  as 


[*n{q}  +  [1>q]T{M  =  {«}  +  {Q(e)} 


(66) 


One  approach  to  the  solution  of  our  system  is  to  form  the  second 
time  derivative  of  Eq. (64) ; 


(67) 


where 


and 


£<h>l 

-{qP(t)}  ] 


aX 

at2 


-fldiXi  +  dij ) 


d  r3h  i 
dt  LdqJ 

[0] 


(68) 
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The  elements  of  the  upper  partition  of  the  matrix  in  Eq.  (68)  are  found 
to  be 


— . .  =  2{  sinCOj-Oi)  +  ±L  -  fj  cos  (0^  —  63^)3  8^ 

dt  or^ 

+  sintOj-ej  +  -  ri  cos  (Bj -6^]  5^  } 


6_dK 
dt  30, 


2[  (  r-rj  +rAXj  )  sinCBj-ej  +  r^  (Gj  —  0± )  cos  (0_j  —  0±)  ]  (  8^  -  8^  ) 


d  3hk 
dt  3z, 


2(zj-zi)(8j,-8i,) 


Writing  Eqs. (66)  and  (67)  together  we  have 


’[*]  [<frq]T" 

r««n  _ 

'  {*}  +  {Qte)} 

[<I>q]  [0] 

M " 

Equations  (69)  can  be  solved  simultaneously  for  the  acceleration 
variables  and  the  Lagrange  multipliers,  to  be  followed  by  numerical 
integration  for  the  generalized  coordinates. 

Constraint  Forces 

To  calculate  the  tether  tension  forces  we  first  establish  the 
corresponding  generalized  force  expressions.  As  was  the  case  with  the 
extensible  tethers,  we  denote  the  force  acting  on  PA  arising  from  its 
connection  to  Pj  by  the  vector  and  let  denote  the  scalar 

tension.  Then  we  can  write 


For  a  tether  in  tension  fi:j>0,  while  fi:j<0  indicates  compression.  If  P± 
and  Pj  are  not  connected,  fi:j=0.  Unlike  the  situation  for  extensible 
tethers,  where  was  expressible  as  an  explicit  function  of  the  system 
state  variables,  here  the  terms  fAj  are  unknowns.  The  total  force  on  Pi 
from  all  its  tether  connections  is,  again, 
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j=l 


The  above  expression  is  identical  to  that  obtained  for  extensible 
tethers.  The  virtual  work  performed  by  the  tether  forces  is  again  given 
by  Eq. (19) ,  which  leads  to  the  generalized  forces  given  in  Eq.  (24)  .  Here 
we  rewrite  Eq. (24)  as 


f.  > 

ril 

(f*  ' 

t  -  t  Eh  Ei2 

-  EiN  ] 

fi2 

=  ir], 

^i2 

,  fiN  J 

where 


Introducing  the  vector  of  tether  tensions;  { T }T  —  T2  Tnt  }  ,  we  can 

write 


*-i2 


v  fiN  J 


=  [B]i {T} 


(71) 


where  [  B  3 ±  is  a  Boolean  selection  matrix  that  establishes  the 
correspondence  between  particle  pairs  and  their  associated  tethers. 
Substituting  Eq. (71)  into  (70)  we  obtain 


Qf  =  inA  tB  ]±  {T}  =  ( A]i  { T } 


The  vector  of  generalized  forces  associated  with  the  tether  tensions  for 
the  entire  system  is  assembled  as 


'  [A],  ‘ 

{Q(1)}  = 

= 

[A]a 

/ - 

IO 

25  2 

_ 

.[A],  . 

{T}  =  [A]  {T} 


(72) 
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In  addition  to  the  tether  tension  forces,  we  consider  external 
constraint  forces  to  act  on  each  particle  to  ensure  satisfaction  of  the 
prescribed  motion  constraints.  These  forces  may  be  thought  to  arise  from 
thrusters.  As  was  shown  in  the  development  leading  to  Eq.  (34)  ,  and 
continuing  the  notation  established  there,  the  generalized  forces 
associated  with  external  constraint  forces  can  be  expressed  as 

{ f(c1  }T  =  (f£>  rlF£>  F<?  j  F<|>  f£  j  -  )  (73) 

It  is  known  from  analytical  mechanics  theory  that  the  generalized 
forces  associated  with  the  constraints  of  Eq. (64)  are  given  by 

{ Q(c)  }  =  -[*q]W  (74) 

The  generalized  constraint  forces  of  Eq. (74)  arise  from  the  inextensible 
tethers  and  the  external  particle  constraint  forces  needed  to  enforce 
any  additional  prescribed  motion  conditions.  In  terms  of  Eqs . (72) - (74) 
we  write 

-[«>q]T{M  =  [A]{T}  +  {F<c)}  (75) 

We  are  interested  in  finding  the  tether  tension  forces  that  most 
closely  satisfy  the  left  hand  side  of  Eq.  (75),  i.e.  we  wish  to  solve  the 
overdetermi*ned  system 

[A]  {T}  =  -[<l>q]T{M 

The  Moore-Penrose  generalized  inverse151  provides  us  with  the  desired 
solution 

{T}  =  -[A  ]+[*q]T{M  (76) 

Substituting  Eq. (76)  into  (75)  yields  a  system  of  particle  forces 
sufficient  to  complete  satisfaction  of  the  constraints 

{F(c)}  =  ([A]  [ A]+  —  [1]  )  [<*>q]T{M 


31 


References 


[1]  Dusto,  A.  R.,  MAn  Analytical  Method  for  Predicting  the  Stability  and 
Control  Characteristics  of  Large  Elastic  Airplanes  at  Subsonic  and 
Supersonic  Speeds",  34th  Meeting  of  the  First  Mechanics  Panel  of  AGARD, 
Marseilles,  France,  April  1969. 

[2]  Canavin,  J.R.,  and  P.W.  Likins,  "Floating  Reference  Frames  for 
Flexible  Spacecraft",  Journal  of  Spacecraft  and  Rockets,  Vol.14,  No. 12, 
December  1977,  pp724-732. 

[3]  Hughes,  P.  C.,  "Spacecraft  Attitude  Dynamics",  J.  Wiley  &  Sons, 

1986,  Appendix  B. 

[4]  Bert,  C.  W. ,  "Material  Damping:  An  Introductory  Review  of 
Mathematical  Models,  Measures  and  Experimental  Techniques",  Journal  of 
Sound  and  Vibration,  Vol .  29,  Number  2,  1973,  pp  129-153. 

[5]  Noble,  B.,  and  J.W.  Daniel,  "Applied  Linear  Algebra",  2nd  Edition, 
Prentice-Hall,  Inc.,  1977,  Chapter  9. 


32 


